Transport in suspended graphene 
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Motivated by recent experiments on suspended graphene showing carrier mobilities as high as 
200, 000 cm 2 /Vs, we theoretically calculate transport properties assuming Coulomb impurities as 
the dominant scattering mechanism. We argue that the substrate-free experiments done in the 
diffusive regime are consistent with our theory and verify many of our earlier predictions including 
(i) removal of the substrate will increase mobility since most of the charged impurities are in the 
substrate, (ii) the minimum conductivity is not universal, but depends on impurity concentration 
with cleaner samples having a higher minimum conductivity. We further argue that experiments on 
suspended graphene put strong constraints on the two parameters involved in our theory, namely, 
the charged impurity concentration nj mp and d, the typical distance of a charged impurity from the 
graphene sheet. The recent experiments on suspended graphene indicate a residual impurity density 
of 1 — 2 x 10 10 cm -2 which are presumably stuck to the graphene interface, compared to impurity 
densities of ~ 10 12 cm -2 for graphene on SiCh substrate. Transport experiments can therefore be 
used as a spectroscopic tool to identify the properties of the remaining impurities in suspended 
graphene. 

PACS numbers: 81.05.Uw; 72.10.-d, 73.40.-c 
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I. INTRODUCTION 



The recent experimental realization of a single layer 
of carbon atoms arranged in a honeycomb lattice has 
prompted much excitement in both theoretical and ex- 
perimental physics communities (For a recent review, see 
Ref. [l] and references therein) . Most experiments so far 
have been performed on a SiC>2 substrate using the me- 
chanical exfoliation technique developed in Ref. [2] fol- 
lowed by optical detection that is sensitive to the thick- 
ness of the substrate. In these experiments one expects 
there to be charged impurities located in the substrate 



with a typical density of n- Y 



10 12 cm 2 with the 



impurities sitting near the substrate-graphene interface 
with a typical distance d < 1 nm from the graphene sheet. 
In previous work^il we argued that removing the sub- 
strate would remove the charged impurities and enhance 
mobility. Other models of scattering such as ripples, on 
the other hand, would lead to enhanced scattering in sus- 
pended graphene thereby reducing the mobility, since a 
freely suspended graphene sheet is expected to be more 
rippled than one glued to a substrate. 

Very recently experiments have been performed on sus- 
pended graphene^ that find remarkably high mobili- 
ties /i = 200, 000 cm 2 /V^s for carrier transport confirm- 
ing that charged impurities in the substrate were indeed 
the dominant scattering mechanism limiting the mobil- 
ity. Signatures from these new experiments (such as the 
conductivity still being linear in carrier density) imply 
that the remaining impurities in suspended graphene are 
long-range and are probably charged dopants lying (i.e. 
stuck) on the graphene surface, or remnants of chemi- 
cals used during the preparation. The minimal model we 
propose here contains only two parameters: the charged 
impurity density and the typical distance of the impurity 



from the graphene sheet. Theoretical estimates for these 
numbers in suspended graphene are hard to come by, for 
example, a recent DFT study-£ argues that potassium and 
sodium dopants are typically located at d ~ 0.3 nm from 
a graphite surface. In what follows we treat "imp and d as 
parameters, and argue that the experimental signatures 
both at high density and low density effectively constrain 
these numbers. Turning the problem around, the recent 
transport experiments on suspended graphene could be 
used as a spectroscopic tool to identify the properties 
of the remaining impurities in suspended graphene. For 
example, a mobility /i = 200, 000 cm 2 /V s, implies that 
"■imp = 1-73 x 10 10 cm -2 , or equivalently, that the dis- 
tance between the remaining impurities is around 76 nm 
with the charged impurities sitting typically at an effec- 
tive distance of 0.1 — 1 nm from the graphene sheet. 



II. THEORETICAL MODEL 



To calculate the transport properties of graphene, we 
consider impurities with an effective 2D density 7ii mp lo- 
cated in a plane at a distance d, either above or below the 
graphene sheet. While in reality, different impurities may 
be located at different distances, one should therefore 
think of d as an effective or typical impurity distance and 
"imp as an average 2D impurity density. Alternatively, 
one could imagine several impurity planes at different 
distances di, and then adding up their respective contri- 
butions to conductivity using Matthiessen's rule. As we 
show below, impurities located farther away are exponen- 
tially suppressed, so in this respect, one could then think 
of d then as distance of the closest charged impurities and 
thus dominating the transport properties. In any case, 
our minimal model for charged impurities contains only 
two parameters 7ii m p and d. 
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For an impurity located at the origin, the bare 
Coulomb potential at the point r in the graphcnc sheet 
is 



Mr) 



Ze 2 



1 



K y/r 2 + d 2 : 



(1) 



where Ze is the impurity charge and k is an effective 
background lattice dielectric constant. All the prop- 
erties of the substrate (or vacuum in the case of sus- 
pended graphene) can be absorbed into a third param- 
eter r s = Ze 2 /nj, where j/H = v-p w c/300. Here c 
is the speed of light, and vp is the constant graphene 
Fermi velocity that relates energy to momentum in the 
linear Dirac spectrum ef = The carrier density 

n = k F /ir. Although we keep r s as an arbitrary number 
in the theory, we note for concreteness that r s — 0.8 for 
graphene on a S1O2 substrate, and r s = 2 for suspended 
graphene (see Refs. HUH for further details on the model). 
The Fourier transform of Eq. [1] gives 



<f>o(q) = 2-K^r s 



4>{q) = 2irjr s 



-qd 



e -qd 



(2) 



where in the second line, we assume linear screening with 



Ak F r s /q 
nr s /2 



q < 2k F 
q > 2fc F . 



(3) 



Although our choice of dielectric function has an artificial 
discontinuity at q — 2fcp, its advantage is that it allows 
for analytic results, and we have checked by numerical 
integration that this choice of dielectric function and the 
Random Phase Approximation^ give identical results 
for all graphene transport properties (see Ref. [H and ref- 
erences therein for a discussion of this point). Further, 
we observe that the chiral properties induced by the sub- 
lattice structure of the honeycomb lattice imply that in 
contrast to most other electronic systems, 2fcp scattering 
is absent in graphene. 

Graphene conductivity can then be calculated in Boltz- 
mann transport theory where a = 2(e 2 / 'h)e-pT ew and 
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|k-k'|)| 2 (l-k.k')|Ckk|^(e fc -£ fc 0, 



where (k is the spinor component of the electron (or hole) 
eigenvector of the 2x2 Dirac Hamiltonian. These expres- 
sions can be simplified by introducing two dimensionless 
variables rj — q/2kp and a — 4kpd — Air l l 2 d^/n to find 



I(a,r s ) 



e 

Ti 

2ri 



I(a,r s 



-ar\ 



drm 2 s/l-r) 2 - , 
o (r? + 2r s ) 2 



(4) 



To calculate the transport properties of graphene from 
Coulomb impurities the integral I(a,r s ) was solved nu- 
merically by AndcS, Cheianov and Falko-P and Hwang 



et ali We explore several limits in which I (a, r s ) can be 
solved analytically. In the "Complete Screening" limit, 
d — » 0,r s — > co, Nomura and MacDonald 11 - showed that 
I(a,r s — > oo) = 7r/32. This result does not capture the 
effect of changing the substrate which in turn changes the 
parameter r s . This complete screening approximation 
is therefore inadequate to study transport in suspended 
graphene even at a qualitative level. For example, as 
we show below, keeping n; mp the same, but just chang- 
ing the dielectric constant from SiC>2 to vacuum, would 
actually decrease mobility by almost a factor of 2. There- 
fore, in suspended graphene the two effects of reducing 
both dielectric constant and n- lmp actually compete to 
determine graphene mobility (see also Ref. fl2l) . Taking 
d — > 0, but keeping the r s dependence, we showed^ that 
7(0, r s ) = 2/G(x = 2r s ) 



G[x] 
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3a; 



3irx 2 x(3x 2 — 2) arccos[l/a;] 
~^2~ + VP = I = 



(5) 



We note here that the following additional limits can be 
solved: I(a,r s — * 0) = (r 2 7r/a)(/i(a) + Li(a)), where I\ 
is a modified Bcsscl function and L\ is a modified Struve 
function. We also have I(a — > 00, r s ) ~ a -3 , where this 
last result says that for very large density, Coulomb im- 
purities in graphene give a(n 3> 1) ~ n 2 ' 5 . However, 
this is not the limit relevant to experiments, which are 
all in the low density regime. The following expansion 
captures both the substrate (i.e. r s dependence) and the 
d dependence, and is relevant for current experiments 
both on suspended graphcnc and for graphcnc on a Si02 
substrate 



a 



h G[x] - dy/nF[x]' 



(6) 



where 

2y/lTX 2 



-4x 2 + 2ttx 3 



x 2 (3 — Ax 2 ) arccos[l/a;] 

v^ 2- T 



(7) 

We note that for most relevant substrates x > 1, (recall 
that for a Si02 substrate x — 2r s w 1 .6 and for suspended 
graphene x ~ 4), and that the ratio A[x] — F[x]/G[x] 
does not change significantly from its large r s asymptote 
of A — > 128/(15^/7?). Therefore for practical compar- 
isons with experiments, one can approximate graphene 
conductivity using 



1 



G[2r s ]m mp J 1-dV^i^ 



(8) 



For Si0 2 substrate G[2 x 0.8 



0.1, giving the 



result of a 



20(e 2 /h)(n/n- m 



and for 



suspended graphene, G[2 x 2] w 0.144 giving a w 
13.86(e 2 //i)(n/rii m p). We find Eq. [5] to be a very use- 
ful formula for all existing bulk graphcnc transport data 
in the literature. 

We now turn to the discussion of mobility /1. Defining 
mobility in graphene is both important and problematic. 
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The importance stems mostly from it being the figure of 
merit for graphene's application as a practical transis- 
tor device. However, close to the Dirac point we have 
V g — > but finite resistivity. If one naively assumes 
that n ~ V g in this regime as is often done, then one 
can get arbitrarily high values of mobility which is com- 
pletely meaningless (a detailed discussion of this point 
can be found in Ref. Il3h . Calculating mobility at very 
high density is complicated by the fact that both super- 
linear contributions (see Eq. [5]for finite d) and sub- linear 
contributions from short-range scatterers (see Ref. HQ] 
and discussion in Sect. IIVI below) come into play, so that 
mobility is not constant for a given sample. In this case, 
mobility would be a non-monotonic function of carrier 
density whose maximum value is given by a competition 
of different scattering mechanisms present in graphene. 
We have consistently argued that the best definition of 
graphene mobility is the slope da/dn just beyond the 
Dirac point conductivity plateau, i.e. only within the 
window of density where the conductivity is strictly lin- 
ear in Vg. We advocate this choice because at low density, 
one can show theoretically that both the sub-linear and 
super-linear contributions to the conductivity are small. 
For this definition of mobility, one can immediately see 
that 



/i[m 2 /y S ] = 



no 



G[x] n 



(9) 



imp 



. Although our theory has two free 
(, lmp and d, knowing the mobility imme- 



where no = 10 10 cm" 
parameters, i.e. n- n 

diately determines ni mp through Eq. [9j A careful char- 
acterization of the super-linear contribution using Eq. [5] 
could then be used to determine d. In addition, as we 
describe below, the low density transport properties i.e. 
the value of the minimum conductivity, the width of the 
minimum conductivity plateau and the offset of the Dirac 
point depend only on these two parameters that in prin- 
ciple could already be determined from the high-density 
transport data. 

So far we have formulated everything in terms of car- 
rier density n — fc|/V. For large external gate volt- 
age, we can map our results to applied gate voltage 
using n = aV g , where a is experimentally determined 
to be 7.2 x 10 10 cm~ 2 V~ 1 for graphene on SiC>2 and 
3.76 x 10 10 cm _2 V _1 for suspended graphene 5,14 . The 
last step is to determine the carrier density n close to the 
Dirac point. The low density conductivity arises from the 
carriers in the charge-impurity induced puddles of elec- 
trons and holesi We argued in Ref. H that at low density 
it was appropriate to use the rms density caused by the 
inhomogeneous potential induced by the same charged 
impurities that are responsible for the high-density trans- 
port properties. A self-consistent calculation gives- 



-1-1 

2 { ' (4 



16E-, \a 



(1 + x)e ax (E 1 [ax] 



(10) 
E 1 [a(l + x)}) 



where Ei(z) — J z °° t^ 1 exp - ' ' dt is the exponential inte- 
gral function and the self-consistency is obtained by set- 
ting a = 4d^/7rn rms . The offset in the Dirac point is 
given by n = n? m ^4n rms (a derivation of this result can 
be found in Ref. y). The value of n rms calculated from 
this self-consistent calculation agrees with a recent cal- 
culation that minimizes the graphene energy functional 
in the presence of disorder— and for d ~ 1 nm is con- 
sistent with the experimental observations of density in- 
homogeneities on a SiC>2 substrate^. Given the choice 
of substrate (which determines r s ) and picking n- lmp and 
d completely determines the conductivity at high den- 
sity and the non-universal disorder dependent minimum 
conductivity plateau at low density. 



III. RESULTS 

To illustrate our results, in what follows, we choose 
the experimental parameters reported in Ref. [H We use 
mobility [i = 30, 000 cm 2 /V^s for graphene on a sub- 
strate (corresponding to n- lm p w 1.7 x 10 11 cm~ 2 ) and 



200,000 cm 2 /Vs for suspended graphene (n 



1.7 x 10 1U cm 



imp 



We also employ a hard cut-off between 



if 



the low and high density regime such that n 
aV g < n lms and n — aV g for aV g > n rms . Shown in Fig.[T] 
is the resistivity as a function of gate voltage for both sus- 
pended graphene and graphene on a SiC>2 substrate. An 
important observation is that the higher mobility sam- 
ples have a lower resistivity which is an important predic- 
tion of our self-consistent theory^. This has been verified 
both in transport measurements on SiC>2 substrate 1 ^ and 
by intentionally adding charged impurities in ultra-high 
vacuum 17 . The recent data 5 on ultra- high mobility sus- 
pended graphene samples further confirm that charged 
impurities provide the dominant scattering mechanism in 
graphene and are responsible for the non-universal min- 
imum conductivity. We note that in plotting resistivity, 
all the super-linear behavior is obscured and the high- 
density curves of the same mobility but different d look 
quite similar. By contrast, from the properties of the low 
density "Dirac plateau" one can extract the magnitude of 
d, but we caution that we only expect our self-consistent 
carrier density theory to agree within a factor of two with 
the experimentally observed plateau width, Dirac point 
offset and the maximum resistivity. A careful analysis 
of the experimental data should be able to extract the 
best values for ni mp and d as has been done previously 
for graphene on a substrat e 14 ' 17 . 

In Fig. [2] we show conductivity as a function of "carrier 
density" defined through n = aVg. For experimentally 
reasonable carrier densities, one can notice that d = 1 nm 
and d = 0.1 nm should be distinguishable by the super- 
linear contributions of the former. In many ways plot- 
ting the data as conductivity instead of resistivity reveals 
much more information. From these results we argue 
that transport data not only provide information about 
^imp but also d, thereby being a tool to not only iden- 
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FIG. 1: (Color online) Transport properties of suspended 
graphene (solid lines) compared to graphene on a SiC>2 sub- 
strate (dashed lines). The self-consistent carrier density the- 
ory predicts that the cleaner suspended graphene samples 
would have a higher minimum conductivity (i.e. lower maxi- 
mum resistivity), smaller off-sets and narrower plateaus. As 
described in the text, we use iHmpiSiOz) = 1.7 X 10 11 cm -2 
(i.e. mobility /i = 30,000 cm 2 /Vs) for graphene on a sub- 
strate and n imp = 1.7 x 10 10 cm -2 ( (i = 200,000 cm 2 /Vs) 
for substrate-free graphene. 
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FIG. 2: (Color online) Suspended graphene conductivity us- 
ing the same values of ni mp as in Fig. Q] and d = 0.1,0.5, 
and 1 nm. Plotting the conductivity instead of resistivity 
reveals some additional properties including the predicted 
super-linear behavior shown in the figure that can be used 
to determine the distance of the charged impurities from the 
graphene sheet. 

tify different types of scattering mechanisms (e.g. short- 
range, long-range, phonons) but for Coulomb impurities, 
it provides information about the location of the impu- 
rities. In Fig. [3] we show our analytical expression Eq. [5] 
(solid lines) compared to a numerical evaluation of Eq. [4] 
(squares) for different values of d and for both suspended 
graphene and on a substrate. One finds that for most of 
the experimental regime, the analytical expression is an 
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FIG. 3: (Color online) Comparison of analytical result Eq. [6] 
(solid line) and a numerical evaluation of Eq. [4] for realistic 
parameters in suspended graphene and on a substrate. The 
curves from top to bottom are suspended graphene with d = 
1 nm and d — 0.1 nm, and graphene on a Si02 substrate with 
d = 1 nm and d = 0.1 nm respectively. The values of 7ii mp 
are the same as in Fig. [1] 



excellent approximation for the Boltzmann conductivity, 
but that at very high-density and large ri, it slightly over- 
estimates the super-linear behavior. However, in this 
high-density regime, one also expects the contribution of 
non-Coulombic short-range scatterers which reduce the 
conductivity. For large d and a particular ratio of short- 
range and Coulomb scatterers, one may even have non- 
monotonic mobility that first increases because of our 
predicted super-linear contribution and then decreases as 
short-range scattering becomes more dominant at higher 
densities. 

Finally, in Fig. [4j we show our calculated carrier mean 
free path (£) as a function of carrier density, comparing 
graphene on Si02 substrate with suspended graphene. It 
is remarkable that, for reasonable parameters (same as 
Fig.[T]) used in these results, the suspended graphene £ > 
1.5 /im for n~ 3 x 10 11 cm -2 whereas for the same den- 
sity range in graphene on Si02 substrate, the charged im- 
purities give rise to an almost constant I ~ 0.1 — 0.2 /im. 
As noted by Bolotin et al£ such long (> 1/im) mean free 
paths are comparable to the 2D system size of the sam- 
ples, implying ballistic transport and a breakdown of our 
diffusive theory. It is important to realize that for mean 
free paths comparable to the sample size, strong short- 
range scattering will be imposed by the sample edges 
leading to a saturation of the conductivity and the mean- 
free-paths around this density. Note, however, that the 
current graphene samples on the substrate are always in 
the diffusive limit, and graphene transport around the 
Dirac point (both with and without substrate) is also 
always in the diffusive limit since the mean free path is 
much shorter than the sample size. Of course, by making 
nanoribbon a 18 ' 19 with sample sizes comparable to £, one 
can artificially induce ballistic transport. 
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FIG. 4: (Color online) Mean free path I for suspended 
graphene (solid lines) and on a substrate (dashed lines) using 
the same parameters as in Fig. [T] For sample sizes L ~ 1/im, 
we observe that close to the Dirac point, transport for both 
suspended graphene and on a substrate is in the diffusive 
regime, where our theory applies. For larger carrier densi- 
ties, suspended graphene crosses over to a ballistic transport 
regime, whereas for graphene on a substrate the transport is 
always diffusive. 



IV. DISCUSSION 

In this work, we have generalized the transport theory 
of Ref. y to the case of suspended graphene^ and de- 
rived analytical results for the super-linear contribution 
to conductivity expected for Coulomb scatterers at a fi- 
nite distance d from the graphene sheet. However, this 
effect becomes pronounced only at high density where 
short-range scattering begins to dominate over Coulomb 
scattering. We note that in the experiments of Bolotin et. 
ali, they were only able to reach rather low carrier den- 
sities, and therefore we believe that Coulomb scattering 
dominates over all other scattering mechanisms in their 
results. However, as was pointed out in Ref.0, there may 
be strong short-range scattering by sample edges when 
the mean free path becomes comparable to the system 
size. However one can anticipate for the Dirac plateau 
in cleaner samples (such as future suspended graphene 
experiments), charged impurities will induce a smaller 
carrier density n rros , even as the ratio n rms /ni mp ~ a 
is increasing. For example, a sample with mobility of 
/i = 500, 000 cm 2 /V~s, would have ra imp = 0.7 x 10 10 cm" 2 
which for d = 0.3 nm gives n rms w 1.2 x 10 10 cm -2 and 
(train ~ 20e 2 /h. This gives I w 600 nm. So despite the 
high-mobility, the peculiar scaling £ ~ kp for Coulomb 
scatterers in graphene implies that close to the charge 
neutral point, we are in the diffusive Boltzmann trans- 
port regime. However, at higher carrier density we have 
£ w 22 kp/riimp, suggesting that the ballistic regime could 



be explored away from the Dirac point. For mesoscopic 
samples constructed in the geometry of £ > L, where L 
is the sample size, our theory does not apply. 

Smaller sample sizes would also increase the effect of 
edge scattering, which act as short-range scattering cen- 
ters. As was argued in Hwang et al4^ at high density 
these begin to dominate showing a non-universal cross- 
over from the low-density regime where a is linear in n 
to the high-density limit where a is constant. So while 
in this work we have focused on Coulomb scattering that 
dominates at low density, probing a wider density range 
could also distinguish between different scattering mech- 
anisms. We emphasize here that detailed measurements 
and analysis of a(n) could thus provide in-depth spec- 
troscopic information about the impurity distribution in 
graphene, far surpassing the resolution of currently avail- 
able scanning microscopy measurements^ which directly 
probe the potential landscape. 

Finally we note that alternative theories for the dom- 
inant scattering mechanism in graphene make predic- 
tions that are inconsistent with the results of Bolotin 
et al£. For example, if ripples were the dominant scat- 
tering mechanism, then suspended graphene would have 
stronger ripples and hence lower mobility. It was one 
of our early predictions that removing charged impuri- 
ties in the substrate would increase mobility and increase 
the minimum conductivity. The initial conventional wis- 
dom in the theoretical community that graphene ex- 
periments were observing the universal Dirac minimum 
conductivity should now be discarded after recent ex- 
perimental developments including intentionally doping 
graphene with potassium in Ref. 17! and measurements 
done on graphene with mobility varying close to three or- 
ders of magnitude in Refs. [5IIT4 and most particularly the 
suspended graphene measurements, all of which clearly 
demonstrate the importance of charged impurity disorder 
controlling the minimum conductivity plateau. It was in 
Refs. that we predicted that charged impurities in 
graphene would give rise to a non-universal minimum 
conductivity whose value would increase with cleaner 
samples. These predictions have now been vindicated in 
light of recent experiments. It, therefore, appears that 
the existing graphene transport experiments arc com- 
pletely dominated by charged impurity scattering, and 
no universal "Dirac point physics" is in play near the 
charge neutrality point. The physics of the charge neu- 
trality point is dominated by diffusive transport through 
electron and hole puddles induced by the charged impu- 
rities. We believe that the close agreement between our 
theory and the transport data in suspended graphene ex- 
periments establishes the dominance of charged impurity 
scattering in graphene beyond any reasonable doubt. 

This work is supported by U.S. ONR and NSF-NRI. 
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